Magnetic trapping of ultracold Rydberg atoms in low angular momentum states 
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We theoretically investigate the quantum properties of nS, nP, and nD Rydberg atoms in a 
magnetic loffe-Pritchard trap. In particular, it is demonstrated that the two-body character of Ry- 
dberg atoms significantly alters the trapping properties opposed to point-like particles with identical 
magnetic moment. Approximate analytical expressions describing the resulting Rydberg trapping 
potentials are derived and their validity is confirmed for experimentally relevant field strengths by 
comparisons to numerical solutions of the underlying Schrodinger equation. In addition to the elec- 
tronic properties, the center of mass dynamics of trapped Rydberg atoms is studied. In particular, we 
analyze the influence of a short-time Rydberg excitation, as required by certain quantum-information 
protocols, on the center of mass dynamics of trapped ground state atoms. A corresponding heating 
rate is derived and the implications for the purity of the density matrix of an encoded qubit are 
investigated. 
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I. INTRODUCTION 

During the past decade, powerful cooling techniques 
enabled remarkable experiments with ultracold atomic 
gases revealing a plethora of intriguing phenomena. 
Among the many fascinating systems are Rydberg atoms 
possessing extraordinary properties [T]. Because of the 
large displacement of the valence electron and the ionic 
core, they are highly polarizable and, therefore, expe- 
rience a strong dipole-dipole interaction amongst each 
other. In ultracold gases, the latter has been shown the- 
oretically [21 [3] and experimentally [H [5l |6l [7l |8] to entail 
a blockade mechanism thereby effectuating a collective 
excitation process of Rydberg atoms [21 UHl E] • More- 
over, two recent experiments demonstrated the block- 
ade between two single atoms a few /im apart [12l [13]. 
The dipole-dipole interaction renders Rydberg atoms also 
promising candidates for the implementation of protocols 
realizing two-qubit quantum gates [2 [5]. A prerequisite 
for the latter is, however, the availability of suitable envi- 
ronments enabling the controlled manipulation of single 
Rydberg atoms and preventing the dephasing of ground 
and Rydberg state. 

Several works have focused on the issue of trapping 
Rydberg atoms based on electric [131 [H], optical [IB] , 
or strong magnetic fields [T71 [H]. Being omnipresent 
in experiments dealing with ultracold atoms, inhomoge- 
neous magnetic fields seem predestined for trapping Ry- 
dberg atoms (even a two-dimensional permanent mag- 
netic lattice of loffe-Pritchard microtraps for ultracold 
atoms has been realized experimentally [IHIEO]). Similar 
to ground state atoms, the magnetic trapping of Ryd- 
berg atoms originates from the interaction of its mag- 
netic moment with the magnetic field. In particular, 
this allows utilizing trap geometries which are well-known 
from ground state atoms. In this spirit, theoretical stud- 



ies recently demonstrated that Rydberg atoms can be 
tightly confined in a magnetic loffe-Pritchard (IP) trap 
fn\ [55] and that one-dimensional Rydberg gases can be 
created and stabilized by means of an additional elec- 
tric field [23, . However, the trapping mechanism relies 
in these studies on high angular momentum electronic 
states that have not been realized yet in experiments 
with ultracold atoms. In a very recent work, the authors 
expanded the former studies to low angular momentum 
nS 1/2 states and showed that the composite nature of Ry- 
dberg atoms, i.e., the fact that they consist of an outer 
electron far away from a compact ionic core, significantly 
alters their trapping properties opposed to point-like par- 
ticles with the same magnetic moment [53]. Furthermore, 
it has been demonstrated how the specific features of the 
Rydberg trapping potential can be probed by means of 
ground state atoms that are off-resonantly coupled to the 
Rydberg state via a two photon laser transition. In the 
present work, we provide a detailed derivation and dis- 
cussion of the Rydberg energy surfaces presented in Ref. 
[51] . Moreover, the trapping potentials arising for the 
nS, nP, and nD states of *^Rb are explored. As we are 
going to show, they possess a reduced azimuthal symme- 
try and a finite trap depth, which can be a few vibrational 
quanta only or less. Choosing the magnetic field param- 
eters appropriately, on the other hand, trapping can be 
achieved with trap depths in the micro-Kelvin regime. 
Implications for quantum information protocols involv- 
ing magnetically trapped Rydberg atoms are discussed. 

In detail we proceed as follows. Section [IT] contains 
a derivation of our working Hamiltonian for low angu- 
lar momentum Rydberg atoms in a loffe-Pritchard trap 
which is solved by means of a hybrid computational ap- 
proach employing basis-set and discretization techniques. 
Section |III| then introduces reasonable approximations 
which allow us to gain analytical solutions for the sta- 
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tionary Schrodinger equation and hence for the trapping 
potentials. In Sec. IV we analyze the resulting energy 



surfaces which serve as a potential for the center of mass 
motion of the Rydberg atom. The range of validity of our 
analytical approach is discussed. Section |V] is dedicated 
to the cm. dynamics within the adiabatic potential sur- 
faces. The question of how the cm. state of a ground 
state atom is altered due to its short-time excitation to 



a Rydberg state is illuminated in Sec. VI A heating rate 
associated with this process is derived. In Sec. |VII| the 
effect of the same process on the purity of the density ma- 
trix of a qubit which is encoded in the hyperfine states 
of a ground state atom is discussed. 



II. HAMILTONIAN 

Along the lines of Ref. [22] we model the mutual in- 
teraction of the highly excited valence electron and the 
remaining closed-shell ionic core of an alkali Rydberg 
atom by an effective potential which is assumed to de- 
pend only on the distance of the two particles. In our 
previous works [STJ |2H [Ml this potential could be 
considered to be purely Coulombic since solely circular 
states with maximum electronic angular momentum were 
investigated. The low angular momentum states of alkali 
atoms, on the other hand, significantly differ from the 
hydrogenic ones because of the finite size and the elec- 
tronic structure of the ionic core. However, the resulting 
core penetration, scattering, and polarization effects can 
be accounted for by employing a model potential of the 
form 



V{r) = Vi{r) 



Zi{r) OLc 



[1 



(1) 



where ac is the static dipole polarizability of the positive- 
ion core while the radial charge Zi{r) is given by 



Zi(r) = 1 + (z - l)e-''i'^ - r(a3 + 04^)6" 



(2) 



where z is the nuclear charge of the neutral atom and is 
the cutoff radius introduced to truncate the unphysical 
short-range behavior of the polarization potential near 
the origin [2F|- Note that Vi{r) depends on the orbital 
angular momentum I via its parameters, i.e., a; = ai(Z) 
and Tc = Tc{f). The resulting binding energies are related 
to the effective quantum number n* and the quantum 
defect 8 W — — ~ 2(n-3)^ D' unless stated 

otherwise, all quantities are given in atomic units. 

The coupling of the charged particles to the external 
magnetic field is introduced via the minimal coupling. 
Pi ^ Pi — 'ZiA(ri), with i e {e, c} denoting the valence 
electron and the remaining ionic core of a Rydberg atom, 
respectively; qi is the charge of the i-th particle and A(x) 
is the vector potential belonging to the magnetic field 
B(x). Including the coupling of the magnetic moments to 
the external field (/Xg and /x^ originate from the electronic 
and nuclear spins, respectively), our initial Hamiltonian 



in the laboratory frame reads (employing ge = ^1, = 
1) 

iJ=i[p, + A(r,)]2 + A_[p^_A(r,)]2 

- • B(re) - /X, • B(r,) + Vi(r) + F,o(L„ S) 

(3) 

with r = |re — rd and M being the mass of the ionic core. 
In contrast to our previous studies [HJ [551 [Ml [23 , one 

has to take into account the fine structure of the atomic 
energy levels: For the magnetic field strengths investi- 
gated in this work, the spin-orbit interaction will lead to 
splittings larger than any Zeeman splitting encountered; 
it is given by 

where L,. and S denote the angular momentum and 
spin of the valence electron, respectively. The term 
[1 — o?V[(r) l'2\~'^ has been introduced to regularize the 
nonphysical divergence near the origin |27j . As usual, 
the field-free electronic eigenstates are labeled by the to- 
tal electronic angular momentum J = -I- S. We re- 
mark that the model potential V/(r) has been developed 
ignoring the fine structure; let us therefore briefiy com- 
ment on the accuracy of Eq. Q in reproducing the fine 
structure intervals. For the 40P state of rubidium, our 
approach yields a good quantitative agreement with the 
experimentally determined fine structure splitting, show- 
ing a deviation of less than one percent. This accuracy 
decreases for higher angular momenta; for the 40Z3 state 
nevertheless a qualitative agreement is found. 

The magnetic field configuration of the loffe-Pritchard 
trap is given by a two-dimensional quadrupole field in the 
2^1, a;2-plane together with a perpendicular offset (loffe- 
) field in the xa-direction. It can be created by sev- 
eral means. The "traditional" macroscopic realization 
uses four parallel current carrying loffe bars which gen- 
erate the two-dimensional quadrupole field. Encompass- 
ing Helmholtz coils create the additional constant field 
|28j . More recent implementations are for example the 
QUIC [29] and the clover- leaf configuration [30]. On a 
microscopic scale, the loffe-Pritchard trap has been im- 
plemented on atom chips by a Z-shaped wire [JT] . The IP 
configuration can be parametrized as B(x) = Bc + B;(x) 
with Be = Sea and B;(x) — G \x\e,\ — ^262]. The corre- 
sponding vector potential reads A(x) = Ac(x) -t- A;(x) 
with Ac(x) = ^ [a;ie2 — 2:261] and A;(x) = Gx\X2G-i, 
where B and G are the loffe field strength and the gradi- 
ent, respectively. The quadratic term Bg oc (a::§ — /9^/2)e3 
that usually arises for a IP configuration can be exactly 
zeroed by geometry, which we are considering in the fol- 
lowing. In actual experimental setups, Bg provides a 
weak confinement also in the a:3-direction. Omitting Bg, 
the magnitude of the magnetic field at a certain posi- 
tion X in space is given by |B(x)| — ^ B'^ + G'^p^ , which 
yields a linear asymptote |B(x)| — > Gp for large coordi- 
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nates {p = yjx\ + x\ ^ B/G) and a harmonic behavior 
|B(x)| « B + iG2p2 close to the origin {p <^B/G). 

After introducing relative and cm. coordinates (r 
and R) |38j and employing the unitary transformation 
U = exp [— |(Bc X r) • R], the Hamiltonian describing 
the Rydberg atom becomes 



UHU^ = Ha 



i[L, + 2S] -B,-/!, -B^R + r) 



He ■ B(R) + A,(R + r) • p + iA,(r)2 + . 

(5) 



Here, Ha = pV2 
nian of an alkali atom possessing the energies E^^\j 



He, 



Vi{r) + T4o(Lr,S) is the Hamilto- 
3n( 

= iAi(R + r)2 + iB,.(r xP) + 
U[Vi{r) + Vso{Lr,S)]W - [Vi{r) + Vso{U,S)] are small 
corrections which can be neglected because of the fol- 
lowing reasons: In the parameter regime we are focusing 
on, the diamagnetic contribution of the gradient field, 
Ai(R + r)^, is small compared to the one of the constant 
loffe field, Ac(r)^. The second contribution of i/corr is 
negligible within our adiabatic approach since (P/M) be- 
comes negligible for ultracold temperatures compared to 
the relative motion (p/m). Finally, the remaining terms 
couple to remote electronic states only and are therefore 
irrelevant. 

The magnetic moments of the particles are connected 
to the electronic spin S and the nuclear spin S accord- 
ing to fi^ ~ — S and /x^ = — ^^S, with gpf being the 
nuclear ^-factor; because of the farge nuclear mass, the 
term involving fi^ is neglected in the following. We re- 
mark that the .Z-component of the cm. momentum com- 
mutes with the Hamiltonian ([S]); hence the longitudinal 
motion can be integrated out by employing plane waves 
\Kz) = exp{—iKzZ). In order to solve the remain- 
ing coupled Schrodinger equation, we employ a Born- 
Oppenheimer separation of the cm. motion and the elec- 
tronic degrees of freedom by projecting Eq. (|5| on the 
electronic eigenfunctions ?/'K(r;R) that parametrically 
depend on the cm. coordinates. We are thereby led to a 
set of decoupled differential equations governing the adi- 
abatic cm. motion within the individual two-dimensional 
energy surfaces i?K(R), i-e., the surfaces i?K(R) serve as 
potentials for the cm. motion of the atom. The non- 
adiabatic (off-diagonal) coupling terms AT that arise 
within this procedure in the kinetic energy term can be 
neglected in our parameter regime since they are sup- 
pressed by the splitting between adjacent energy surfaces 
[21j . As will be shown in Sec. Ill the latter is propor- 
tional to the loffe field strength B, i.e., the non-adiabatic 
couplings are proportional to powers of l/B. 

The electronic eigenfunctions and energies are found by 
a standard basis set method utilizing the field-free eigen- 
functions \k) — \njmjls) of Ha whose spin and angular 
parts [jirijls) are given by the spin-orbit coupled gener- 
alized spherical harmonics yj^mj ,i |32j . For the radial de- 
gree of freedom, a discrete variable representation (DVR) 
based on generalized Laguerre polynomials is employed 




FIG. 1: (Color online) Schematic overview of the energy scales 
and coupling terms involved. Starting from the degenerate 
hydrogen energy spectrum on the left, it is shown how the 
quantum defect - modeled by the potential Vi{r) - separates 
the low angular momentum states. The spin-orbit coupling 
Vao(Lr,S) then yields the fine structure splitting for fixed I. 
Within a given j-manifold, the Hamiltonian Ha + |(?jmj |B| 
resembles the coupling of a point-like particle to the mag- 
netic field B(R). The two-body character of the Rydberg 
atom, which is represented by H' , only contributes if energet- 
ically remote levels are considered as well: it admixes states 
of different n, I, j, and nij thereby qualitatively changing the 
shape of the surfaces. 



|33j . The latter provides a non-uniform grid for the ra- 
dial coordinate which is more dense close to the origin 
and hence especially suited for representing radial Ryd- 
berg wave functions. Since in the DVR scheme the poten- 
tial matrix element evaluation is equivalent to a Gaussian 
quadrature rule, representing the Hamiltonian (|5| - espe- 
cially Vi{r) and the derivative terms arising from the mo- 
mentum operator p - becomes particularly efficient. The 
numerical diagonalization of the resulting Hamiltonian 
matrix (in the limit P — > 0) then yields the electronic 
eigenfunctions -0/1 (r;R) and energies i?K(R) which both 
parametrically depend on the cm. coordinates R. Con- 
vergence is ensured by appropriately choosing the size of 
the field-free basis as well as the underlying DVR grid 
size. 



III. ANALYTICAL APPROACH 

While the above described numerical treatment of the 
electronic Hamiltonian offers accurate results, we derive 
in this section analytical but approximative expressions 
for the adiabatic energy surfaces, which provides us with 
a profound understanding of the underlying physics. We 
start by considering only a single fine structure man- 
ifold, i.e., fixed total angular momentum j for given 
I. Such an assumption is motivated by the fact that 
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TABLE I: Coefficients of the linear fit of Ci{n) = cf +Cf 
in the range 35 < n < 45 for the nS, nP, and nD states of the 
*^Rb atom. Note that negative magnetic quantum numbers 
nij yield the same results as their positive counterparts and 
are consequently omitted. The fitted d are calculated using 
Eq. with n' G [n - 10, n + 10]. 





State 






rW 






'S'l/2 


TUj = 


1/2 


-0.4813 


-0.00027 


-0.4813 


-0.00027 


Pl/2 


mj = 


1/2 


-0.4484 


-0.00148 


-0.4484 


-0.00148 


P3/2 


nij = 


1/2 


-0.4541 


-0.00152 


-0.4316 


-0.00164 


P3/2 


ruj = 


3/2 


-0.4391 


-0.00160 


-0.4616 


-0.00149 


D3/2 


TUj = 


1/2 


-0.4570 


-0.00069 


-0.4326 


-0.00011 


D3/2 


ruj = 


3/2 


-0.4407 


-0.00030 


-0.4652 


-0.00088 


D5/2 


ruj = 


1/2 


-0.4570 


-0.00073 


-0.4287 


-0.00006 


D5/2 


ruj = 


3/2 


-0.4500 


-0.00057 


-0.4429 


-0.00040 


D5/2 


ruj = 


5/2 


-0.4358 


-0.00023 


-0.4712 


-0.00107 



the fine structure dominates over the Zeeman splitting 
for the field strengths we are interested in, cf. Fig. [T| 
Within this regime, the contribution A;(R + r) • p = 
G{XY + Xy + xY + xy)pz of Hamiltonian ([5| can be 
simplified as follows. We rewrite 



1 z t 

yPz = ^{yPz- zpy)- -[yz, Ha] + -[yz,Vi{r) + Vso{'Lr,S)] 

. (6) 

bearing in mind that only the action of any involved 
operator within a single j-manifold is considered. The 
first commutator in the above equation then vanishes 
due to the degeneracy of the eigenenergies of Ha and 
since no coupling to different n-, 1-, or J-states is consid- 
ered. Likewise, the second commutator in Eq. ([6| van- 
ishes since neither Vi{r) nor VsoiJ^r,^) depend on the 
magnetic quantum number rrij. Hence, within a given 
j-manifold, the electronic part of Hamiltonian (|5| can be 
approximated by 



He^HA + - [Lr + 2S] • B(R) + GXYp, + Hr , (7) 



where we substituted ypz 



1 T 

2^x 



and similarly xpz 



-\Ly. The contribution Hr = A;(r) • p + B((r) • S + 
|Ac(r)^ only depends on the relative coordinate and - 
as we will show later - for a wide range of field strengths 
can approximately be regarded as a mere energy offset to 
the electronic energy surfaces; we will restrict ourselves 
to this regime and hence omit Hr in the following. 

The first two terms of Hamiltonian ^ can be diago- 
nalized analytically by applying the spatially dependent 
transformation 



(8) 



-GX 



tion; 7 and /? denote the rotation angles: 

-GY 

sin7 : 
sin P - 
cos 7 : 
cos (3 -■ 



+ G2X2 ' 
^/B^ + G^X^ 

B 



(9) 
(10) 

(11) 
(12) 



+ ■ 

The transformed Hamiltonian becomes 

,J,^B^ + G^{X^ + Y^),+H' (13) 



UrH^Ul = Ha 



2 3j •- 



with gj ~ 



s{s+l)-l[l + l) 



, UrHA0 = Ha, and H' 



2 ^ 2i(j + l) 

GXYUrPzUj. . Like for ground state atoms, the second 
term of Eq. ( 13 ) represents the coupling of a point- like 



particle to the magnetic field via its magnetic moment 

As depicted in Fig.[T| H' couples to different n, I, j, and 
rrij and hence vanishes within one j-manifold. The first 
two terms of Eq. ( |13| , on the other hand, are diagonal, 
giving rise to the electronic potential energy surface 



2 9 3' 



,VS2 + G2(X2 + r2) (14) 



for a given electronic state \k) — \njmjls). Note that 
such a state refers to the rotated frame of reference. Only 
there, rrij constitutes a good quantum number; in the 
laboratory frame of reference to,- is not conserved. The 



surfaces Eq. (14 1 are rotationally symmetric around the 
Z-axis and confining for rrij > 0. For small radii {p = 
y/X^+Y^ <C B/G) an expansion up to second order 
yields a harmonic potential 



1 



E(^Hp)^Ef + y,m,B+-Mu;'p' 



(15) 



with the trap frequency defined by lo 



G 



we find a linear behavior E)^'\p) fa E'^ + ^gjUijGp when 
the center of mass is far from the Z-axis {p ^ B/G). In 
the harmonic part of the potential, the cm. energies are 
thus given by 



Ej^\0) + {ly + l)uj , 1^ = + uy e N (16) 



with a splitting of lo between adjacent cm. states; see 
Sec. |V] for a more detailed discussion. The separation 
between adjacent electronic energy surfaces at the origin, 
on the other hand, is given by AE^ = \gjB. The size 
of the cm. ground state (y — 0) in such a harmonic 
potential evaluates to [p) = sjir /2^/ Mw [22] . 



that rotates the z-axis into the local magnetic field direc- 



The remaining term H' of Hamiltonian (13 1 can be 
treated perturbatively. While it vanishes in first order. 
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FIG. 2: (Color online) Range of validity for deriving Eq. (24 1, 
(a) 2G^XVB^ < 1 resulting in X < B/V2G. In the figure, 
the value X+/(i3/\/2G) is shown which should be ^ 1. Fig- 
ure (b) shows G^X^jB^ which should be ^ I as well. In both 
cases I = 2,j = rrij — 5/2 is used. 



second-order perturbation theory yields 

£;(2)(R) = G'X'Y' y \(^\UrP^^rW)\' (17) 



(18) 



where k = \njmjls) are the eigenstates of th e trans 



formed Hamiltonian UrHJJl — iJ', cf. Eq. (13 1. Since 



i?i°''(R) resembles the confinement of a ground state 

atom, we attribute i?„ (R) to the composite nature of 
the Rydberg atom, i.e., the fact that it consists of a Ryd- 
berg electron far apart from its ionic core. Like the mag- 

f 21 

netic field itself, E)^ ' (R) shows no continuous azimuthal 
symmetry but rather a discrete one. 



Equation (18) is derived by employing UrPzU} — 
Px sin/? — sin 7 cos /3 -1-^2 cos7COS/3 and p = i[i?A,r] — 
i\yi{r) + Vso,y] ~ i[i?A,r] [Ml- Expanding the modu- 



lus square in Eq. (18 1, one obtains mixed terms of the 
form (kIj/Ik') -|- Employing the 

standard basis of spherical harmonics and consequently 
using {k,\x\k') g M as well as {Hi\y\K') — — (kIj/Ik')*, this 
sum vanishes. The matrix element of z obeys a different 
selection rule, namely. Am; = opposed to Ato/ = ±1 of 
X and y; hence, mixed terms involving (k|z|k') vanish as 
well. Consequently, only the matrix elements |(K|a;|K')p, 
and remain and the second order 

energy contribution can be parametrized as 

+ Cy sin^7 cos^/3 + cos^7 cos^/3) (19) 
Cz 



1 



(sin^/3 -I- sin^7COS^/3) 



(20) 



E 



(21) 



where C^, = Cy = C^, = C for I = and C^, 



Cy oth- 



erwise (since |(«;|a;|/t') | — \{K\y\K')\). Note that the pa- 
rameters Ci depend on the state k under investigation. 
Since E^'' — E'^, oc and |(K|a:i|K')p oc n^, a hnear 

(2) 

scaling of E)i (R) with the quantum number n is antic- 
ipated, i.e., Ci{n) ~ c|°^ + C'j^^^n. Resulting from a fit 
of calculated Ci values within the range 35 < n < 45, in 
Tab. |l] the coefficients C^^'^ are tabulated for the nS, nP, 
and nD states of the ^^Rb atom. All considered states 
show a similar behavior: The magnitude of Ci is close 
to —1/2 and shows a rather weak n-dependence. In par- 
ticular, Cx ' 
cf. Eq. ([Ml 



Cz and therefore E^^\yi) « • G^X^F^ 
We remark that for smaller n, Ci{n) de- 
viates from the linear behavior in favor of a more rapid 
decrease. 

In the last part of this section, let us reconsider the 
adiabatic energy surfaces for the cm. motion, including 
now the contribution of H' . That is, we investigate the 
approximate, but analytical solutions 



i?.(R)^4")(R) + £;(2)(R)-i?: 



(22) 



of Hamiltonian ([7|). In particular, we concentrate on the 

diagonal of the surfaces {X — Y) where E)^ (R) is max- 
imal. The approximation Cx = Cz (which is exact for 
nSi/2 states) then yields 



E^iX = Y) 



i^B^ + 2G^X^ + CzG^X'^ , (23) 



which shows only a local minimum at the origin since the 

(2) 

surface drops off for large cm. coordinates when E^, (R) 
dominates (note that Cz < 0), see also Fig. [3] The po- 
sitions of the maxima which enclose the minimum are 
approximately given by 



X^ 



± 



91 

52 



4.BCz 
9jmj 



-1/2 



±Xn 1 - 



G2 
2S2 



:Xr 



with the length scale Xq 



(24) 

I only depending on 



the loffe field strength. The depth of the potential well 
associated with the minimum correspondingly evaluates 
to 



E^{X 

1 



Y^X, 



G'^Xl 
B 



1 - 



-EM 

C^Xj 
B2 



+CzG^X^ 1-2 



52 



(25) 



(26) 



Note that the first approximation in Eq. (24 1 holds for 
2G^X^/B^ < 1 and the second one for G^X^/B'^ < 
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FIG. 3: (Color online) First Row: Contour plots of the electronic potential energy surfaces -Ek(R) of the 405*1/2 (first column), 
40Pi/2 (second column), 4OP3/2 (third column), 4O-D3/2 (fourth column), and 4OD5/2 (fifth column) states with rrij = j for 
the magnetic field configuration B = 1 G, G = 2.5 T/m. Second Row: Section along X = Y oi the same energy surfaces; 
the contribution of Ei:^\'R) (dashed lines) is shown in addition. Third Row: Depth of the potentials as given by Eq. (251 



as a function of the magnetic field configuration. For all subfigures, the energy scale is given in terms of the trap frequency 
u! = G^gjmj/2MB. 



1. The corresponding range of validity is illustrated in 
Fig.|2] For a lofFe field strength of B = 1 G, the above ap- 
proximations hold for gradients up to 10 T/m; at higher 
B even larger gradients are eligible. 



IV. TRAPPING POTENTIALS 

In the following section we are going to discuss the cal- 
culated electronic potential energy surfaces for the nS, 
nP, and nD states of the *^Rb atom in detail. In par- 
ticular, the range of validity of the above derived an- 
alytic expression [Eq. (22l] is demonstrated. As a gen- 
eral example, we address the magnetic field configuration 
B — 1 G, G = 2.5 Tm~^, which yields a trap frequency 
of a; = 27r X 319 Hz. A similar field configuration is also 
found in current experiments |34j . 

In Figure |3] the electronic potential energy surfaces 
E,{R.) of the 405i/2, 40Pi/2, 4OP3/2, 407^3/2, and AOD^/^ 
states with rrij = j are illustrated. In addition, also sec- 
tions along X = Y of these surfaces are provided. On a 
first glance, the energy surfaces originating from different 
electronic states seem to differ quite substantially. How- 
ever, qualitatively they are very similar, as we are going 
to argue in the following. For all surfaces presented in 
Fig. |3] the contribution of the composite character of the 



metry of Ej^^R) into a four-fold one. Moreover, the in- 
terplay between the harmonic confinement E^\'R) and 

the unbounded contribution E^ (R) gives rise to a fi- 
nite trap depth along the diagonal X = Y; see second 
row of Fig. 3 Since the coefficient C, of eI?\r) IS ap- 
proximately of the same magnitude for all states consid- 
ered, cf. Tab. |T] the trap depth depends on the magni- 
tude of the magnetic moment /j c>c gjJ. Consequently, 



the j 



= 1 + 1/2 electronic states show a deeper 



confinement than their j 



1 — 1/2 counterparts 



(2) 

Rydberg atom, i.e., i?K (R) changes the azimuthal sym- 



and the depth increases further with increasing orbital 
angular momentum I. For the examples given in Fig. [3] 
this means that the quadratic approximation to the trap- 
ping potential for the 405*1/2 state is already violated at 
about two oscillator energies, while for the 40Z?5/2 it is 
fine up to 10 w. This trend is confirmed in the third row 
of Fig. |3]where the depth of the potential as a function of 
the field configuration is displayed: For the 4OD5/2 state, 
the trap depth easily exceeds 100 cu within the given pa- 
rameter range, while in the case of the 405*1/2 state there 
is a substantial regime of field strengths where not a sin- 
gle center of mass state can be confined, i.e., the trap 
depth being <liu. Nevertheless, also for the 405*1/2 Ry- 
dberg state the field parameters B and G can be adjusted 
such that trapping is possible, i.e., the trap depth being 
much larger than the trap frequency. Similarly, the trap- 
ping potential of the 4OD5/2 Rydberg state can be chosen 



Kl -3 -2 -lX[/im]l 2 3 
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FIG. 4: (Color online) Sections along X — Y of the energy 
surfaces of the (a) 4051/2, (b) 4OP3/2, and (c) 4O-D5/2 states 
for the field configuration B — 0.1 G, G = 10 T/m, which 
yields a trap frequency of u/y/^gjmj — 2-k x 4 kHz. 



very shallow by going to sufficiently strong lofFe fields. 
We remark that the results presented i n Fig. [3| are g iven 
in units of the trap frequency u> = Gyjgjrrij /2MB; the 
latter holds, of course, only near the origin. For larger 

radii, the contribution Ek. (R) flattens the potential re- 
sulting in smaller trap frequencies and hence in a higher 
number of center of mass states that can be confined. 
Moreover, for very high gradients the harmonic expan- 
sion of the magnetic field strength becomes progressively 
worse. 

As can be deduced from the third row of Fig.|3] increas- 
ing the relative strength of the field gradient, i.e., either 
increasing G directly or decreasing the offset field B for 
fixed G, leads to a larger number of bound center of mass 
states ~ independently of the state under consideration. 
However, since the anti-trapping contribution E^, (R) 
quadratically increases with the field gradient G, we ex- 
pect this trend to reverse for sufficiently high gradients. 
Indeed, for a loffe field oi B — 1 G the trap depth starts to 
decrease for field gradients G > 200 Tm"^ for B = 0.1 G 
this trend already starts at G > 5Tm"^. Similarly, for a 
fixed field gradient together with a decreasing offset field 
B we find a decrease of the trap depth if B < 0.15G or 
B < 0.03 G for G = 10 Tm"^ and G = lTm"\ respec- 
tively. 

In the following, let us investigate the question if elec- 
tronic energy surfaces belonging to different states inter- 
sect each other. In Refs. |2TJ [52], where high angular 
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FIG. 5: (Color online) (a) Section along X = F of the energy 
surface of the 40-D5/2,mj —5/2 state for the field configura- 
tion B = 1 G, G = 100 T/m which yields a trap frequency of 
u = 2n X 22.1 kHz near the origin. The solid line is yield by 
the numerical diagonalization of Hamiltonian (|5| , the dashed 
one by Eq. (22 1. The surfaces have been offset to zero at the 



origin, (b) Difference (in terms of the trap frequency to) be- 



tween the analytic expression Eq. ( 22 \ and the result of the 
numerical diagonalization of Hami Ionian ([sjl . (c) Same as 
(b) but for y = rather than X = Y. 



momentum states are considered, this issue is essential: 
there, the high level of degeneracy of the system leads 
to non-adiabatic crossings of the surfaces. As a conse- 
quence, only the circular electronic state {nii = I = n— 1) 
provides stable trapping. For the low angular momentum 
states of ^'^Rb we are considering here, however, the fine 
structure splitting for different j and the varying quan- 
tum defect for different I separate the energy surfaces by 
lifting the degeneracy, therefore preventing their cross- 
ing. As a result, it is sufficient in our case to investigate 
the energy surface spectrum for fixed j, i.e., only as a 
function of the magnetic quantum number. In Figure [4] 
sections along the diagonal of the energy surfaces of the 
multiplets of the 405*1/2, 4OP3/2, and 40Z35/2 states are 
presented. In order to show a strong spatial dependence, 
we choose an extreme case concerning the ratio of the 
loffc field compared to the gradient field, namely, B — 0.1 
G, G=10T/m. Even for such a high gradient (of course 
much higher gradients can be achieved on atom chips), 
the energy surfaces remain well separated with a mini- 
mum distance of \gjB at the trap center. Hence, each 
surface can be considered separately for trapping and our 
adiabatic approach is not limited by non-adiabatic inter- 
actions. Note that the anti-trapping of raj < states is 

even enhanced by the contribution (R)- 

The above investigations employ the analytical expres- 
sion Eq. ( 22 1 rather than the exact numerical solutions 
of Hamiltonian Hence an estimation of the range 
of validity of our results is necessary. To this end, we 
provide in Fig. [sja) a comparison between the analyti- 
cal expression according to Eq. (|22[ ) and the numerical 
diagonalization of Hamiltonian (|5| for the 'extreme' field 
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configuration of i? = 1 G, G = 100 T/m (the quanti- 
tative agreement improves for smaller gradients / larger 
lofFe fields). As one can observe, even for such a strong 
gradient Eq. (22 1 yields satisfactory results: The devia- 



tion within the spatial range considered is less than 0.2 ui, 
cf. Fig.jsjb), iij being the splitting of the cm. states and 
representing the smallest energy scale of the system. We 
remark that Fig. shows results along the diagonals 
X — Y of the surfaces where the deviation is at maxi- 
mum. However, even along the axes, where Er{R) van- 
ishes, a perfect agreement cannot be found, cf. Fig.jsj^c). 
This residual deviation is due to the purely electronic 
terms Hr of Hamiltonian ^ which have been neglected 
in deriving Eq. ( 22 1 . Although Hr does not explicitly de- 



pend on the cm. coordinate R, it introduces an implicit 
cm. dependency by changing the electronic state. This 
can be easily understood in the rotated frame of refer- 
ence, i.e., after applying the unitary transformation Ur'. 
There one has to consider UrHrUl which introduces a 
cm. dependency exphcitly. We stress that for gradients 
weaker and/or loffe fields stronger than in Fig. [s] Ur is 
closer to unity and therefore the contribution of Hr be- 
comes even less important in these cases and hence can 
be neglected. 



V. CM. WAVE FUNCTIONS 

In this section we discuss the eigenfunctions x(R-) of 
the cm. Hamiltonian 



Hr 



2M 



+ Sk(R) 



(27) 



where -Ek(R-) are the previously calculated potentials. In 
particular, we are going to elucidate the differences to 
the harmonic oscillator eigenstates, that are yielded by 
considering solely the 'unperturbed' potential Ej^\'R) 
(thr ough out this section, the harmonic approximation 

Eq. (151 for the potential £^i°''(R) is assumed). The en- 
ergies and eigenstates are then computed using second 
order perturbation theory. We remark that the validity 
of the harmonic approximation together with the use of 
perturbation theory has been ensured by comparing with 
the results obtained by the numerical diagonalization of 



Hamiltonian (27 1 



Before presenting our results, let us comment on the is- 
sue of the finite radiative lifetime of Rydberg atoms which 
might spoil the experimental observation of the cm. mo- 
tion. The lifetime can be parametrized as r = r'(n — (5^)^ 
where one finds r' — 1.43 ns and 7 = 2.94 for I = 0, 
t' = 2.76 ns and 7 = 3.02 for Z = 1, and r' = 2.09 ns 
and 7 = 2.85 for ? = 2 [3F. For the 405i/2 Rydberg 
state, this yields a radiative lifetime of r = 58 /is. If we 
compare this to the typical time scale — 2tt/uj of the 
cm. motion, one finds that for the envisaged field con- 
figuration (B — 1 G, G — 2.5 T/m) r^^ = 3 ms is orders 
of magnitudes larger than the radiative lifetime, which 
renders the resolution of the cm. motion experimentally 



impossible. This drawback can be alleviated by several 
means. First of all, one can consider higher principal 
quantum numbers n which increases the lifetime substan- 
tially. For example, the 605*1/2 state already possesses a 
radiative lifetime of r = 206 fis. The changes in the trap- 
ping potential, on the other hand, are marginal as can be 
seen from the weak n-dependence of the coefficient , cf . 
Tab.|l] note that E^\'R,) is n- independent. Additionally 
to increasing n, one can augment the trap frequency by 
increasing the gradient field and/or decreasing the loffe 
field (which might necessitate atom chip traps [5T1 [55] V 
As an example, the field configuration B — 0.1 G and 
G = 50 T/m yields — 50 fis. Furthermore, one might 
also employ the Rydberg states which possess a 

longer lifetime (r = 155 /xs and r = 0.5 ms for n = 40 
and n = 60, respectively) and at the same time cause 
a higher trap frequency {t^ = 2 ms and = 35 /is for 
5 = 1 G, G 2.5 T/m and S = 0.1 G, G = 50 T/m, 
respectively) . 

As an illustrative example, let us investigate again the 
405*1/2 Rydberg state combined with the magnetic field 
parameters i? = 1 G and G = 2.5 T/m in the following, 
despite the above mentioned restrictions. In this case, 
the resulting trapping potential -Ek(R) confines only a 
very limited number of cm. states, namely, twelve. Con- 
sequently, already low cm. excitations show an apprecia- 
ble deviation from the harmonic behavior which makes 
the influence of the perturbative effects of E)^ ' (R) par- 
ticularly visible. For the case of £;i°^(R) and small cm. 
radii one yields a harmonic potential. In this case, the 
Hamiltonian (27 1 decouples in X and F, i.e., the to- 



tal cm. wave function can be written as a product of 
two independent harmonic oscillator states in X and Y: 
X(R) = X^xi'JR-) = X^A^) ■ Xuy{y)- As a consequence, 
the corresponding energies only depend on the sum of 
the individual cm. excitations v = 
{v + l)-fold degeneracy. 

In anticipation of considering (R) as well, it is 
advisable to employ adapted eigenstates which account 
for the symmetry of Hamiltonian (27) including 



Vy and show a 



(2) 

Ek (R)- The decomposition of such symmetry adapted 
eigenstates in terms of the product states Xu^^y^) can 
be found in the fourth column of Tab. |lT] together with 
their corresponding symmetry label given in the second 
column. Note that these states are still degenerate in 
case of the potential e'"^\r), cf. third column of Tab. II 

The inclusion of E^. (R) lifts this degeneracy by mix- 
ing states of equal symmetry according to the vanish- 

f2l 

ing integral rule |37|: since E^ (R) is of Ai symme- 
try, i.e., being totally symmetric the cm. matrix element 
(2) 

{x'\Ek Ix) is only non- vanishing if |x) and |x') possess 
the same symmetry. Moreover, the perturbation of the 
form ^ X^Y^ yields the selection rules Av^ e {0, ±2} 
and Avy e {0,±2}. Inspecting the wave functions as 



obtained by diagonalizing Hamiltonian ( 27 1 within a for 



merly degenerate i^-manifold, both the symmetry con- 
straints as well as the selection rules become apparent; 
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see sixth column of Tab. |TT1 

The energies i?^™ of the first 12 eigenstates are tabu- 
lated in Tab.|ll]for both potentials Ej^\R) ("h.o.", third 
column) and eI^\r) + E^^\K) ("perturbed", fifth col- 
umn). While eI^\'R) yields energies E^ = {v+ 1)lu with 
^' = 0, 1, 2, . . . (we assumed a perfectly harmonic poten- 
tial), the eigenenergies belonging to Ej^\'R) + Ej^\'R) 
deviate from this rule: as the harmonic potential is flat- 
tened by the contribution E^ (R) , the energies are below 
the harmonic ones. The remaining degeneracies which 
appear for odd j/ (e.g., states 6-9) can be explained by 
the symmetry properties of the involved states: in this 
case, only E symmetry is encountered. The latter has 
a two-dimensional irreducible representation hence the 
appearance of degenerate pairs. 

Finally, let us briefly comment on the issue of tunnel- 
ing. States that are conflned within the potentials shown 



in Sec. IV may escape the trap by tunneling through the 
potential barrier along the diagonals. While this pro- 
cess most certainly plays no role for configurations where 
the time scale of the trap frequency is large compared 
to the radiative lifetime, a priori it is not clear if tun- 
neling becomes crucial for tighter traps. For this rea- 
son, we estimated the lifetime associated with the tun- 
neling process by investigating the transmission proba- 
bility for the ith excited cm. state in one dimension, 

Pt = exp(-2 ^2M[E,{X = Y) - E^^"^ dX), where 

the integration limits a and b are determined by the con- 
dition E^{X = Y) = E'^"^. Since the Rydberg atom 'hits' 
the potential barriers twice per trapping period, the loss 
rate can be roughly estimated by 2u}Pt . Actual values of 
Pt for the cm. states discussed in this section are given 
in the last column of Tab. For tighter magnetic traps, 
where more cm. states can be confined, Pt substantially 
decreases further. Hence, tunneling only has to be con- 
sidered for highly excited cm. excitations close to the top 
of the barrier. 



VI. PARAMETRIC HEATING 

Utilizing state-dependent (Rydberg-Rydberg) interac- 
tions for quantum information protocols necessitates the 
excitation of trapped ground state atoms to a Rydberg 
state by a 7r-pulse [Sj |3j . When the excitation process is 
much shorter than the timescale of the external motion, 
such an excitation effectively causes a sudden change 
of the trapping potential. This couples and thus redis- 
tributes the initial cm. quantum state to neighboring 
levels which, in general, increases the cm. energy (hence 
we will denote this process as "parametric heating" in 
the following). In this section, we investigate this effect 
and calculate the corresponding heating rates. 

Suppose we have a ^''Rb atom in its 5S'i/2, F — mp — 2 
electronic ground state which is at t = instantaneously 
excited to the Rydberg state 405*1/2, ^Tij = 1/2 and af- 



ter a short period of time t' again de-excited to its elec- 
tronic ground state. Furthermore, we assume the atom 
to reside in a well defined cm. state at < = 0, i.e., 
x(R,i = 0) = Xv^yy{^)\ note that Xv^Vy{^) denote the 
cm. cigenfunctions of the ground state atom rather than 
the Rydberg atom. Except for the contribution En (R), 
both electronic states give rise to the same trapping po- 
tential Ef\li) HO], i.e., in the simplest approximation 

[which is neglecting E)^ (R)] the cm. state is not affected 
by the excitation to the Rydberg level. If we account 

(2) 

for the extra term E^, (R), on the other hand, the sit- 
uation changes substantially. We consider the sequence 
ground state — > Rydberg state — > ground state, where all 
transitions are carried out by fast 7r-pulses. E^^(R) can 
then be considered as a perturbation of the ground state 
trapping potential which acts for the time interval dur- 
ing which the atom resides in the Rydberg level, i.e., 
Q < t < t' . As shown in Sec. |v] E^^\ll) mixes cm. 
states according to the selection rules ^i^x/y — and 
^v^jy = ±2; hence the Rydberg excitation leads for 
the ground state atom to the admixture of lower- and 
higher-lying cm. levels with v' = v, v' = v :k2, and 

= 1/ ± 4, where v = Vy. Note that we adopt here 
again the approximation of a purely harmonic potential 
E^^\r) cx \Mu?{X'' + Y^). 

Within time-dependent perturbation theory, the prob- 
ability of a transition \VxV,^ Wx^'y) of the cm. state 
of a ground state atom due to its short-time Rydberg 
excitation is consequently given by 



with f{t',uj) 



'M\E'^\R)\ux^y) f{t',Q) (28) 



sin^(cit72) 



and & ^ {ly' — ly) LO — Av ■ uj 

[32] . The average rate to make a transition to state \v'^v'y) 
within the time interval t' consequently reads 



R,j 



(29) 



This allows us to define a heating rate as 



(30) 



{ 



sm^(2u)t') 



4^2 



yy + 1) {vx - 



where we used the recurrence relation 

2 En:iWvx{Vx - l)5 ^,^^^^2 + {"^.Vx H 

^J{vx + '^){,Vx + 2)5v'^y^+2] of the 

(2), 



)+ 



1, 3 
2^+4 



(31) 



harmonic oscilla- 



tor cigenfunctions and assumed Ek "^(R) = C^G X Y . 
Note that E^,^y^ > independent of the initial state, i.e., 
cooling is not possible. For short times t' <C 1/w, one 
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Hamiltonian (27 1 for the harmonic potential Ek'^TV) 



Fourth Column: Corre- 



TABLE II: Third Column: Energies of the cm. 
sponding symmetry adapted eigenstates. Fifth Column: Second order perturbation theory energies including the contribution 
eI^\'R). Sixth Column: Corresponding zero order perturbation theory wave functions. The symmetries of the states encoun- 
tered are listed in the second column. Seventh column: Tunneling probability Pt for the energies given in the fifth column. 



state no. 


symmetry 


h.o. 

energies 


symmetry adapted 
eigenstates 


perturbed 
energies 


zero order 
eigenstates 


tunneling 
probability 


ground 




UJ 


Xoo 


0.9859a; 


Xoo 


3.19 X 10"** 


1 


\ J 




Xoi 


1.9566 a; 


Xoi 


1.35 X 10"^ 


2 




Xio 


1.9566 a; 


Xio 


1.35 X lO"" 


3 


B2 




Xii 


2.8665 a; 


Xii 


3.38 X 10~^ 


4 


Ai 




(X02 + X2o)/\/2 


2.8942 a; 


(X02 + X2o)/\/2 


3.71 X 10"^ 


5 


Bi 




(X02 - X2o)/V2 


2.9572 a; 


(X02 - X2o)/V2 


4.60 X 10^^ 


6 


{4 


Au 


X12 


3.7490 a; 


0.937x12 + 0.349x30 


6.44 X 10"^ 


7 


4uj 


X21 


3.7490 a; 


0.937x21 + 0.349x03 


6.44 X 10"^ 


8 


{-} 


4uj 


X03 


3.9156 a; 


0.937x30 - 0.349x12 


0.11 


9 


4uj 


X30 


3.9156 a; 


0.937x03 - 0.349x21 


0.11 


10 


Ai 


5 U! 


X22 


4.5616a; 


0.937x22 + 0.257(x4o + X04) 


0.87 


11 


B2 


5 U! 


(Xl3 + X3l)/V2 


4.5684 a; 


(Xl3 +X3l)/^ 


0.89 


12 


A2 


5 U! 


(Xi3 - X3i)/\/2 


> AE^ 


(Xi3 - X3i)/\/2 




12 


Bi 


5 U! 


(X04 - X4o)/\/2 


> AE^ 


(X04 - X40)/^ 




12 


Ai 


5 ui 


(X04 + X4o)/y2 


> AE^ 


0.663(xo4 + X4o) - 0.349x22 





can approximate 
linear increase of the heating rate in time. 



In Figure [oj the parametric heating E^^^^t' in terms of 
the trap frequency uj and as a function of the Rydberg 
excitation period t' is illustrated for several cm. initial 
states and magnetic field configurations for the Rydberg 
state 405*1/2. As one can observe, the heating mainly 
depends on the loffe field strength B rather than on 
the magnetic field gradient G. An increase of the latter 
barely changes E^^^jj! juj while a stronger loffe field re- 
sults in a substantial increase. As expected from Eq. (31 1, 
Ey^^^t' jio also significantly increases if the ground state 
atom is initially in an excited cm. state. However, for 
the given examples the overall heating within the radia- 
tive lifetime of the Rydberg atom turns out to be very 
moderate with E^^^Jl < la; (for B —\Q and G — 2.5 
T/m, E'jy^jy t' — lu! corresponds to 15 nK). Hence, only 



t'^ which gives an overall In this case, the heating rate reads 



for high cm. levels ly 



and long times t' the 



above described excitation of the cm. motion of an ultra- 
cold sample of Rb atoms due to the Rydberg excitation 
is expected to become an issue. 



Finally, let us briefly comment on what is expected 
for a thermal atom where the cm. state is not a pure 
state but rather a mixture according to the Boltzmann 
distribution /^(T) = g^e-(''+i)'^/'=''^/Z(r), Z{T) = 
Yl°^=o9v^~~'^'^'^^^^^''^'^ being the partition function and 
gu — V + ^ the degeneracy of the i^th excited cm. state. 



coth^ 



2M\ 
VlClG^ 



2 khT 



''( 



(32) 

(33) 
(34) 



Equation ( 33 ) is obtained by approximating 



sin^(a;i') 



t'"^ for short times t' ^ 1/ui and further simplified to 
Eq. (|34[ ) by assuming k\jT ^ uj. As expected from 
Eq. ( |3l| , E(T) rapidly increases with the temperature 
T since higher cm. excitations are populated. 



VII. DEPHASING 

Besides the parametric heating due to the short-time 
Rydberg excitation of a ground state atom - as discussed 
in the previous section - the dephasing of the cm. motion 
of the Rydberg and the ground state might become an is- 
sue for experimental schemes realizing quantum informa- 
tion protocols. Let us consider the situation as described 
in Rcf. [5], i.e., we have two ground states denoted by 
|0) and |1) where only the latter is coupled to a Rydberg 
state \r) by a laser transition. The density operator of the 
internal degree of freedom, i.e., only considering the elec- 
tronic state, of such a two state system can generally be 
written as Pint =a|0)(0|-H(l-a)|l)(l|-H6|l)(0|-|-6*|0)(l| 
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FIG. 6: (Color online) Parametric heating Ev^^yt' /ui of a 
trapped ground state atom as a function of the time being 
excited to the Rydberg level 40Si/2. Several initial cm. states 
and magnetic field configurations are considered: B = 1 G, 
G = 2.5 T/m (solid), B = 1 G, G = lOT/m (short-dashed), 
and B = lOG, G = 2.5 T/m (dashed-dotted) for u^, = Vy = Q 
as well as B = IG, G = 2.5T/m (dotted) and B = lOG, 
G = 2.5 T/m (long-dashed) for = Vy = 2. 



giving rise to the density matrix 

^'"^^(fol-a)- (3^) 

If we assume furthermore that both ground states are 
identically prepared with respect to their external, i.e., 
cm. motion, the total density matrix factorizes into an 
internal and external contribution, ptot = Pint ® Poxt, 
where pcxt = "^^yVy Vv^vy 

For various implementations of quantum information 
protocols now the Rydberg state |r) comes into play. 
Suppose that state |1) is excited to |r) for a given time 
i! . As pointed out in Sec. 



VI 



this will influence its 
cm. motion by causing transitions \vxV^ —* where 

|X> ^ E.>;^^^S(i')l^>;>; denotes the am- 

plitude for being at time t' in state \v'^v[^ if initially 
residing in state \v.j.v^. Hence, after the short-time Ry- 
dberg excitation of solely state |1), the density matrix 
does not decouple anymore and consequently reads 

t^^t^y 

+ {l-a)\l){l\®\x){x\ 
+ b\l){Q\®\x){yxUy\ 

+ b*\Q){l\®\v.,Vy){x\]. (36) 

Any qubit-related measurement, however, only acts on 
the internal degrees of freedom, i.e., the electronic states. 
As a consequence, the relevant object in this case is the 
reduced density matrix where the cm. degree of freedom 



is traced out. Defining (3 
eventually yields 

TrextPtot = 



Y.>.^u,.P'^.'^y{^xVy\x), one 



Comparing this result to the case where the internal and 
external degree of freedom factorize, Trext(/Oint ® Pext) = 
Pint, it is clear that the short-time Rydberg excitation 
will inevitably influence the properties of our system. 
In order to quantify this effect, we consider the purity 
7'(p) — Trp^ of the reduced density matrix. In particu- 
lar, if the system is initially prepared in a pure internal 
state (as for example (|0) ± |l))/\/2, which is envisaged 
for the realization of a two qubit CNOT gate [12]) the 
above described process is expected to decrease its pu- 
rity. Indeed, one yields 

^(TrextPtot) = a2 (1 - af + 2\b\''\P\'' (38) 
= 7'i„t-2|6|2(l-|/3|2) (39) 

where V-mt = a? + {l — a)'^ + 2\h\^ denotes the purity of the 
system only considering the internal degree of freedom. 
Note that for any pure state Vmt = 1 is found. Hence, the 
reduction of the purity is determined by the magnitude 
of and therefore by the overlap integrals {vxyy\x) of 
the cm. wavefunction after the Rydberg excitation. 

A particularly illustrative situation arises, if the atom 
is initially prepared in its cm. ground state, i.e., Pv^vy — 
5ou^5ouy and |/3p = |(00|x)p correspondingly. In this 
case, |/3p is given by the probability of finding the atom 
still in the cm. ground state after being excited to the 
Rydberg state for the time t' . According to Sec. |VI[ we 
find 



= 1 - 



'2/^4 



4^2 



~ 1 - t''^ 

and therefore 

'P(TrcxtPtot) « "Pint 



2\h\ 



(40) 

(41) 
(42) 

(43) 



a b*(3* 
bp 1-a 



(37) 



As expected, the decrease of the purity depends explicitly 
(for short times t' even quadratically) on the time t' of 
being excited to the Rydberg level. 



VIII. CONCLUSION 

We theoretically investigated the quantum properties 
of Rydberg atoms in a magnetic loffe-Pritchard trap. In 
particular, the electronic properties and the center of 
mass dynamics of the low angular momentum nS, nP, 
and nD states of ^^Rb have been studied. It turns out 
that the composite nature of Rydberg atoms, i.e., the fact 
that it consists of an outer electron far away from a com- 
pact ionic core, significantly alters the coupling of the 
electronic motion to the inhomogeneous magnetic field 
of the loffe-Pritchard trap. We demonstrated that this 
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leads to qualitative changes in the trapping potentials, 
namely, the appearance of an de-confining contribution 
which reduces the azimuthal symmetry to C^y. As a 
consequence, the resulting energy surfaces - which char- 
acterize the trapping potentials - possess a finite depth. 
Analytical expressions describing the surfaces were de- 
rived and the applicability of the applied perturbative 
treatment has been verified for experimentally relevant 
field strengths by comparison with numerical solutions 
of the underlying Schrodinger equation. Exemplary en- 
ergy surfaces of the fully polarized n = 40, I = 0,1,2, 
rrij = j states for the magnetic field configuration B = I 
G, G = 2.5 T/m were provided. A clear deviation from 
the harmonic confinement of a point-like particle with 
a trap depth of only a few vibrational quanta could be 
observed. Choosing difi^erent magnetic field parameters, 
on the other hand, trapping can be achieved with trap 
depths in the micro-Kelvin regime. The non-harmonicity 
of the Rydberg trapping potential becomes also apparent 
in the resulting center of mass dynamics: The additional 
contribution due to the two-body character of the Ryd- 
berg atom mixes the "unperturbed" harmonic eigenstates 
and thereby partially lifts their degeneracy. For an atom 
in its electronic ground state that is excited to a Ryd- 
berg state only for a short period of time, this provides 
a mechanism for parametric heating by populating ex- 
cited center of mass states. The corresponding heating 



rate as a function of the initial center of mass state of the 
ground state atom has been derived. In the framework of 
quantum information protocols involving the short-time 
population of Rydberg atoms, it has been demonstrated 
that the same mechanism can lead to a decrease of the 
purity of the involved qubit states. 

A rather natural extension of the present work would 
be the investigation of magnetic field configurations other 
than the loffe-Pritchard trap. For example, it is expected 
also in the case of a three-dimensional quadrupole field 
that similar terms associated with the composite nature 
of the Rydberg atom arise, significantly altering the trap- 
ping potential compared to the point-like particle de- 
scription. 
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